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ABSTRACT 


An  analytical  model  and  computer  simulation  are  presented  for  a  radome 
located  in  the  near  field  of  an  antenna.  Using  the  computer  code  described  here, 
design  tradeoffs  can  be  performed  between  electrical,  stractural,  and  aerodynamic 
properties  of  the  radome.  The  code  is  based  on  a  method  of  moments  solution  to  the 
E-field  integral  equation  for  bodies  of  arbitrary  shape.  Measured  radiation  patterns 
for  AGM-88  High  Speed  Antiradiation  Missile  (HARM)  and  AIM-9C  missile 
radomes  are  compared  to  computed  data. 
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L  INTRODUCTION 


The  radome  is  an  essential  part  of  communications  and  radar  systems  deployed  on 
aircraft.  It  provides  environmental  protection  for  the  antenna  and  therefore  must  withstand 
extreme  aerod5aiamic  and  thermal  stresses  (such  as  on  a  supersonic  missile)  and  must  also 
effectively  transmit  electromagnetic  radiation.  For  example,  the  radome  is  a  critical 
conponent  on  an  anti-radiation  missile  (ARM).  An  ARM  uses  a  passive  antenna  to  home  in 
on  enemy  radar  emissions.  The  emitted  energy  passes  through  the  missile’s  radome,  which 
has  some  loss  and  also  introduces  beam  pointing  error.  The  main  beam’s  field  of  view  (FOV) 
must  be  narrow  enough  for  proper  target  discrimination,  but  wide  enough  for  off-boresight 
intercept.  The  radome  can  also  increase  side  lobe  levels,  which  must  remain  low  enough  to 
prevent  azimuth  and  elevation  errors.  The  result  is  degraded  antenna  performance  when  the 
radome  is  present  relative  to  when  the  antenna  is  operating  isolated  in  free  space. 

Traditionally,  aerod5Tiamic,  structural,  and  thermal  requirements  have  dictated  the 
design  of  a  radome,  with  the  electrical  performance  taken  as  a  consequence.  The  analytical 
model  and  conputer  code  described  in  this  thesis  can  be  used  to  optimize  the  overall  design 
of  the  radome  since  electrical  design  tradeoffs  can  be  performed  concurrently  with  the 
aerodynamic  and  structural  engineering.  Therefore,  an  accurate  analytical  model  of  the 
antenna  and  radome  is  required  for  efficient,  cost-effective  radome  design. 

There  are  basically  two  methods  of  analyzing  radomes  of  arbitrary  shape:  (1)  ray 
tracing  (geometrical  optics),  and  (2)  the  method  of  moments  (MM).  For  ray  tracing,  a 
spherical  wave  is  assumed  radiated  from  the  antenna,  which  is  approximately  planar  at  the 
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radome  so  that  the  Fresnel  reflection  coefficients  can  be  used.  Rays  are  traced  from  the 
antenna  to  the  radome  surface,  where  reflections  and  refractions  occur.  The  transmitted  rays 
are  used  to  reconstruct  an  aperture  distribution.  The  far  field  is  computed  using  radiation 
integrals  [Ref.  1]. 

Ray  tracing  is  only  applicable  for  electrically  large  radomes  that  are  in  the  far  field  of 
the  antenna.  For  low  side  lobe  antennas,  miiltiple  reflections  are  important,  and  a  significant 
amount  of  bookkeeping  is  required  to  determine  the  ray  trajectories.  An  example  of  ray 
tracing  is  illustrated  in  Figure  1.1. 


Figure  1.1  Ray  tracing 
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The  computer  code  presented  here  uses  a  method  of  moments  solution  to  the  electric 
field  integral  equation  to  calculate  the  effect  of  a  dielectric  radome  shell  in  the  near  field  of 
an  antenna.  The  method  of  moments  is  rigorous  in  that  all  multiple  reflections  between  the 
radome  walls  are  included.  As  will  be  demonstrated,  the  code  is  flexible  as  well  as  accurate 
and  handles  radomes  of  any  shape  and  size  (with  size  limited  only  by  the  available  computer 
memory). 

Chapter  II  reviews  the  analysis  methods  and  computer  codes.  An  overview  of  the 
method  of  moments  is  presented.  Essentially  the  electric  field  integral  equation  is  reduced  to 
a  matrix  problem  which  is  readily  solved  by  numerical  methods.  Also,  the  equivalent  single¬ 
layer  impedance  is  derived  for  multiple  layer  radome  walls.  The  impedance  is  calculated  for 
two  test  radomes:  an  AGM-88  HARM  and  an  AIM-9C.  The  impedance  values  are  used  by 
the  method  of  moments  computer  codes.  Chapter  11  concludes  with  a  discussion  of  these 
codes.  RADOME,  a  code  developed  by  Francis  [Ref.  2]  and  later  improved  by  Klopp 
[Ref  3],  is  limited  to  bodies  of  revolution  (BOR),  and  uses  a  thin-shell  approximation  for  the 
radome  wall  In  this  thesis,  a  code  developed  by  Sandia  National  Laboratories  called  PATCH 
[Ref.  4]  is  applied  to  the  radome  problem.  It  represents  the  radome  as  a  collection  of 
triangular  patches  which  allows  for  virtually  any  geometry  as  well  as  various  dielectric 
materials.  Cartesian  coordinates  describe  antenna  and  radome  geometry  in  PATCH,  while 
spherical  coordinates  describe  the  electromagnetic  field.  The  coordinate  systems  are  shown 
in  Figure  1.2.  Typically  the  radome  is  oriented  with  the  axis  of  symmetry  coincident  with  the 
z-axis. 
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Chapter  HI  explains  the  procedure  used  to  modify  and  validate  PATCH.  An  antenna 
aperture  model  is  developed  for  a  standard  gain  horn  that  is  used  for  both  PATCH  and 
RADOME  programs.  A  test  case  was  run  and  the  computed  results  from  both  programs  are 
conpared  to  measured  data  obtained  in  the  anechoic  chamber.  An  additional  test  case  is  used 
to  demonstrate  the  capability  of  PATCH  to  accurately  calculate  multiple  reflections. 

Chapter  IV  presents  the  data  obtained  by  analyzing  the  two  test  radomes.  The 
anechoic  chamber  facility  is  described  and  measuring  techniques  are  discussed.  Several 
inprovements  have  been  made  since  the  chamber’s  inception.  Automated  measurements  are 
made  using  a  HP-8510  network  analyzer  linked  to  a  personal  computer. 

Finally,  Chapter  V  presents  a  summary  of  results  and  conclusions  based  on  working 
with  the  code  as  a  radome  analysis  tool.  Sources  of  measurement  error  are  discussed  and 
recommendations  are  made  to  enhance  the  existing  computer  models  and  codes. 
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n.  ANALYSIS  METHODS  AND  COMPUTER  CODES 


A.  METHOD  OF  MOMENTS 

The  computer  codes  that  will  be  discussed  in  the  following  sections  use  the  method 
of  moments  (MM)  to  calculate  radome  effects.  A  detailed  discussion  of  this  method  can  be 
found  in  [Ref.  1].  Only  an  overview  of  the  method  is  presented  here. 

The  electric  field  integral  equation  (EFIE)  is  commonly  encountered  in  electro¬ 
magnetics.  This  equation  relates  the  scattered  electric  field  to  the  electric  current  on  the 
surface  of  an  object.  The  EFEE  is  classified  mathematically  as  a  Fredholm  equation  since  the 
limits  of  integration  are  fixed  (i.e.,  over  the  body).  Because  the  unknown  variable  appears 
only  in  the  integrand,  it  is  a  Fredholm  equation  of  the  first  kind. 

To  derive  the  EFIE,  we  first  define  the  coordinate  system  shown  in  Figure  2.1.  The 
target,  or  scattering  body,  is  assumed  (temporarily)  to  be  a  perfect  electric  conductor  (PEC). 
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We  will  introduce  a  correction  term  later  that  will  remove  the  PEC  restriction.  Note  that  in 


general,  S  includes  both  the  radome  and  antenna  surfaces.  We  define 
=  |r  -  =  ^(x  -  x^'  +  (y  - y^'  +  (z  - 


and  the  firee-space  Green’s  fimction  as 


G(F,T') 


e 

4tiR 


(2.1) 


(2.2) 


After  some  manipulation  [Ref  1,  pp.  117-118],  the  final  form  of  the  EFIE  suitable  for 
solution  by  a  numerical  technique  is 


tan 


o)pJ/F)G(F,F0 


(2.3) 


where  £.  is  the  incident  electric  field,  is  the  surface  current  density,  p  is  the  permeability, 
e  is  the  dielectric  constant,  and  o)  is  the  angular  frequency.  Since  the  total  electric  field  is  the 
sum  of  the  incident  and  scattered  electric  fields,  by  imposing  the  PEC  boundary  condition 
(total  tangential  electric  field  is  zero).  Equation  2.3  states  that  the  tangential  component  of  E. 
equals  the  tangential  component  of  the  scattered  field. 

This  form  of  the  equation  can  be  applied  to  any  arbitrary  surface.  The  incident  field 
is  always  known,  and  the  surface  current  density  is  the  only  unknown.  The  incident  field  can 
be  a  plane  wave  or  an  impressed  voltage  on  the  surface  of  the  body.  Potentials  on  the  antenna 
surfaces  produce  an  electric  field  which  then  becomes  the  incident  field  on  a  scattering  body 
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such  as  a  radome.  This  interaction  between  the  antenna  and  radome  is  embedded  in 
Equation  2.3. 

The  method  of  moments  is  a  procedure  used  to  reduce  the  integral  equation  to  a 
matrix  problem.  We  start  by  representing  the  current  in  Equation  2.3  as  a  series  with 
unknown  expansion  coefficients,  /„ . 

h  =  (2.4) 

n=l 


The  vectors  are  known  as  the  basis  functions.  These  are  either  subsectional  or  entire 
domain  functions,  and  are  chosen  to  be  accurate  representations  of  the  current  and 
mathematically  convenient  (easy  to  integrate  and  differentiate).  Typical  subdomain  basis 
functions  include  impulse  functions,  rectangular  pulses,  and  overlapping  triangular  pulses. 
Entire  domain  functions  include  sinusoids  and  exponentials.  A  testing  procedure  is  performed 
in  which  both  sides  of  Equation  2.3  are  midt^lied  by  each  testing  function  and  integrated  over 
its  domain.  Thus  N  equations  are  obtained,  each  with  N  terms.  The  simultaneous  set  of 
equations  can  be  cast  into  matrix  form  and  when  solved  yield 

I  =  Z-iV  (2.5) 

where  I  is  the  vector  containing  the  expansion  coefficients,  /„ .  Z  is  an  impedance  matrix  and 
V  the  excitation  vector.  Generally  both  of  these  have  complex  elements.  The  specific  form 
of  the  elements  depends  on  the  choice  of  basis  functions.  The  coefficients  from  Equation  2.5 
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are  then  used  in  Equation  2.4,  which  yields  the  surface  current,  .  This  in  turn  is  used  in  the 

radiation  integrals  to  find  the  fields  due  to  the  surface  current. 

The  element  of  the  excitation  vector  V  in  Equation  2.5  has  the  form 

V.  - 1 1  E,-y: ‘Is 

where  Jj"  is  the  testing  function  and  E.  the  incident  antenna  field.  S„  denotes  the  portion 
of  the  radome  surface  over  which  is  nonzero.  In  general,  the  can  lie  in  the  near  field 
of  the  antenna  and  therefore  no  far  field  approximations  can  be  applied.  The  antenna  models 
in  PATCH  and  RADOME  differ,  and  both  will  be  described  later.  Neither  code  has 
restrictions  on  the  location  of  the  radome  relative  to  the  antenna. 

Equation  2.3  applies  to  a  perfectly  conducting  surface.  For  thin  dielectric  shells  such 
as  radomes,  the  EFDE  for  perfect  conductors  is  modified  by  adding  a  load  term  outside  of  the 
brackets  [Ref  1,  p.  158].  This  thin-shell  approximation  replaces  the  volume  polarization 
currents  flowing  throughout  the  dielectric  with  a  thin  sheet  of  surface  current.  The  result  is 
that  the  elements  of  Z  in  Equation  2.5  are  also  modified  by  the  addition  of  a  load  term. 

Overall,  the  method  of  moments  provides  a  rigorous  solution  for  the  induced  current 
density  on  a  body.  When  properly  applied,  the  series  representation  for  the  current  will 
converge  to  the  actual  current  as  the  number  of  basis  functions  is  increased  [Ref  1].  Since 
the  solution  reduces  to  a  matrix  problem,  MM  is  well-suited  for  computer  coding. 

Two  MM  solutions  have  been  applied  to  the  radome  problem  in  this  thesis.  One 
solution  (RADOME)  applies  only  to  bodies  of  revolution.  By  exploiting  the  rotational 
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symmetry  of  bodies  such  as  ogives  and  hemispheres,  basis  functions  can  be  chosen  that 
reduce  the  number  of  linear  equations  that  must  be  solved.  The  second  solution  (PATCH) 
uses  triangular  subdomains,  which  are  the  most  flexible  for  modeling  arbitrary  surfaces. 

B.  RADOME  WALL  MODEL  AND  IMPEDANCE  CALCULATION 

The  impedance  of  a  multi-layered  (sandwich)  radome  can  be  calculated  using  a 
stepped  transmission  line  analogy.  The  impedance  of  the  region  on  the  exterior  of  the 


r-n 

Zin 


Figure  2.2  Stepped  impedance  transformer 

radome,  ,  is  transformed  through  each  layer  until  the  region  on  the  interior  side  is  reached. 
This  transformed  inpedance,  ,  is  the  impedance  that  is  presented  to  a  plane  wave  incident 
on  the  radome.  The  schematic  representation  of  a  sandwich  radome  is  shown  in  Figure  2.2. 
To  apply  the  transmission  line  model  to  the  radome,  we  first  define  the  quantities  shown  in 
Figure  2.3.  Each  layer  n  has  a  thickness,  ,  and  a  relative  dielectric  constant,  q,  .  The 
exterior  region  has  the  highest  index  and  the  interior  region  the  lowest.  The  exterior  is  region 
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number  N,  so  there  are  N-2  layers  of  radome  material.  Let  r|„  be  the  intrinsic  impedance  of 
the  layer  n  and  Z„  be  the  wave  impedance  looking  into  the  layer  n  at  an  incidence  angle  0„. 
For  perpendicular  polarization 


COS0 

n 


(2.7) 


and  for  parallel  polarization 


=  'n„cos0„ . 

ft  ‘ft  ft 


(2.8) 


is  the  transformed  iirqjedance  that  is  seen  looking  through  a  sandwich  of  layers  as  shown 
in  Figure  2.3.  The  transformed  impedance  is 


^  ^ntanh(Y„g 
^  2„,itanh(Y„t„). 


1  ^n^(iV-l) 


(2.9) 


where  as  in  [Ref.  5] 


(2.10) 


(2.11) 


€ 


e'  -  je'' 


(2.12) 
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Figure  2.3  Multi-layered  radome  wall  geometry 


The  input  impedance  is  computed  by  transforming  the  free  space  load  of  the  exterior 
back  to  the  interior  surface.  An  input  impedance  is  computed  at  each  interface  and  then  used 
as  a  load  to  be  transformed  back  to  the  previous  interface  [Ref.  6].  This  is  repeated  until 
is  computed  for  the  final  layer. 

The  change  in  incidence  angle  at  the  interface  between  materials  with  different 
dielectric  constants  (refraction)  can  be  determined  by  using  a  form  of  Snell’s  Law.  This  is 
expressed  as 


sin0 


n+l 


—  sin6_ 


or 


sin0 


sin0.  . 


(2.13) 
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If  necessary,  an  equivalent  single  layer  radome  (denoted  by  the  subscript  e)  can  be 


derived  from  a  multiple  layer  radome  by  using  the  equation 


+  Z^tanh(Y^t^) 
+  Z„tanh(Y,i,). 


(2.14) 


Since  Z;„  is  complex,  Z^  and  the  product  Ye4  can  be  determined. 

Once  known,  the  impedance  Z^  can  be  used  to  solve  for  the  reflection  coefficient  seen 
by  an  incident  wave 


P  = 


Z.  -  Z 

in  o 

z.  +  z 

in  o 


(2.15) 


When  Z,„  and  ^  are  equal,  p  =  0,  and  all  energy  is  transmitted.  When  there  is  a  large 
difference  between  Z;„  and  ^ ,  |  p  |  ->  1 ,  and  a  large  reflection  results.  Hence,  the  best 
radomes  have  impedances  that  are  close  to  the  intrinsic  impedance  of  free  space, 
approximately  377  ohms. 

The  PATCH  and  RADOME  programs  require  radome  impedance  as  an  input  parameter. 
For  a  sandwich  waU,  the  value  can  be  calculated  from  Equation  2.9.  The  impedance  is  allowed 
to  change  as  a  function  of  position  on  the  radome,  but  it  is  assumed  constant  with  angle.  Using 
a  corrputer  program  (see  Appendix  C)  to  calculate  Z^  for  arange  of  incidence  angles,  an  estimate 
of  the  radome  impedance  was  obtained  that  was  relatively  constant  over  a  broad  range  of 
incidence  angles.  An  average  value  of  impedance  was  chosen  for  each  of  the  radomes  measured 
and  used  in  the  computer  simulations. 
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Two  particular  radomes  were  used  for  testing  and  analysis  in  this  thesis:  (1)  an 
AGM-88  HARM,  and  (2)  an  AIM-9C.  The  cross  sectional  geometry  for  the  walls  are  shown 
in  Figures  2.4  and  2.5,  respectively.  The  computer  generated  impedances  for  the  HARM 
radome  are  shown  in  Figure  2.6.  The  impedance  is  relatively  constant  over  the  range 
0°  <  0  <  50® .  From  the  geometry  of  the  HARM  radome  and  the  antenna,  it  was  determined 


Figure  2.5  AJM-9C  radome  wall 
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Figure  2.6  HARM  radome  impedance  as  a  function  of  incidence  angle 


that  the  average  incidence  angle  along  the  surface  of  the  radome  was  approximately  35°. 
Thus,  an  impedance  of  278  -  y’Sl  ohms  was  used  in  the  computer  simulations.  For  the 
AIM-9C  radome,  the  impedance  calculation  was  simplified  due  to  the  presence  of  just  one 
layer  and  nearly  normal  incidence  around  the  hemispherical  surface.  Using  Equation  2.9,  the 
AIM-9C  impedance  was  234  -  yi66  ohms.  The  irrpedance  values  for  both  radomes  are  listed 
in  Table  2.1. 


Table  2.1  Radome  Impedance  Used  in  the  Computer  Calculations 


Radome 

Material 

Complex  Impedance 

HARM 

Fiberglass  composite 

278  -  7*3 1  ohms 

AIM-9C 

Alumina 

234  -  yi66  ohms 
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C.  PATCH  AND  RELATED  PROGRAMS 


The  program  PATCH  was  produced  by  Sandia  National  Laboratories  in  1988 
[Ref.  4],  It  is  a  FORTRAN  program  that  uses  the  method  of  moments  to  calculate 
electromagnetic  fields  from  scattering  and  radiating  bodies.  PATCH  performs  the 
electromagnetic  calculations,  while  a  collection  of  pre-  and  post-processing  programs  are 
available  for  generating,  formatting,  converting,  and  error  checking  the  input  and  output  of 
the  PATCH  program.  A  brief  description  of  each  program  is  given  in  Table  2.2.  The 
programs  and  their  relationships  are  shown  m  Figure  2.7.  Most  of  the  programs^  were 
generated  at  the  Naval  Postgraduate  School  to  allow  automatic  meshing  of  bodies  built  in 
AutoCAD  and  to  translate  the  meshed  surfaces  to  PATCH  input  format. 

During  execution,  PATCH  reads  an  input  file  with  geometry  and  calculation 
parameters.  Either  of  two  commercial  conputer  aided  drawing  (CAD)  programs,  AutoCAD 
or  ACAD,  can  be  used  to  specify  a  shape.  These  programs  are  extremely  versatile,  and  are 
unlimited  in  the  scope  of  geometries  that  can  be  generated.  Files  can  be  transferred  between 
the  two  using  the  IGES  format.  Whether  the  shape  is  generated  in  AutoCAD  or  ACAD,  the 
later  program  is  used  to  constract  a  shell  mesh,  which  is  written  in  ACAD’s  “facet”  format. 
Meshing  results  in  a  discretized  body  composed  of  triangular  facets.  For  example,  a  simple 
rectangular  plate  is  shown  in  Figure  2.8.  The  degree  of  fineness  (or  coarseness)  of  the  mesh 
is  determined  by  specifying  the  number  of  facets,  or  edges,  along  the  horizontal  and  vertical 


‘Commercially  available  program  names  are  in  ALL  CAPITALS,  locally  generated 
programs  are  in  boldface,  and  inpul/output  files  are  in  italics.  Locally  generated 
programs  have  various  authors,  and  have  been  modified  to  suit  individual  users. 
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Table  2.2  Summary  of  PATCH  and  related  processing  programs 


Name 

Brief  Description 

Source 

AutoCAD 

CAD  drawing  program 

Autodesk 

ACAD 

CAD  program  capable  of  auto  meshing 
surfaces 

Lockheed-Martin 
(Ft.  Worth,  TX) 

pta.x 

PATCH-to-ACAD  data  translator 

NPS 

atp.x 

ACAD-to-PATCH  data  translator 

NPS 

knitx 

Checks  PATCH  data  files  for  duplicate  edges 
and  nodes 

NPS 

bldmatx 

Generates  MATLAB  “.m”  files  for  plotting 
geometry  using  PATCH  data  files 

NPS 

pitpatch.m 

A  MATLAB  program  that  plots  the  geometry 
files  generated  by  bldmatx 

NPS 

combine.m 

A  MATLAB  program  that  combines  separate 
geometry  files  into  one  object  for  use  in 

PATCH 

S.M.  Herzog 

BUILDN5.X 

PATCH  preprocessor  to  generate  input 
geometry  files  and  specify  calculation 
parameters 

Sandia  National 
Laboratories 

PATCH 

(patch2v.x) 

Numerical  electromagnetics  code  that 
computes  scattered  and  radiated  fields  using 
the  method  of  moments. 

Sandia  National 
Laboratories 

sides.  As  a  rale  of  thumb,  the  lengths  of  the  individual  edges  should  be  approximately  one 
tenth  of  a  wavelength  (0.1  A).  For  larger  bodies,  this  means  a  large  number  of  elements,  and 
a  correspondingly  larger  facet  file. 

The  facet  file  generated  by  ACAD  must  then  be  converted  to  PATCH  text  format. 
This  is  accomplished  by  the  program  atp.x  (ACAD  to  PATCH).  The  program  pta.x 
performs  the  reverse  data  translation.  The  PATCH  version  of  a  facet  file  defines  the 
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Figure  2.7  Relationships  between  PATCH  and  its  pre-  and  post-processing  programs 
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Figure  2.8  Plate  composed  of  triangular  facets 


geometry  in  terms  of  nodes  and  edges.  Nodes  are  listed  with  their  respective  Cartesian 
coordinates,  and  edge  segments  are  listed  with  their  two  defining  nodes.  An  example  of  this 
will  be  presented  later.  The  output  of  atp.x  is  used  in  the  program  BUILDN5.X,  which 
constracts  the  final  ASCII  input  file  that  PATCH  will  read.  When  BUILDN5.X  is  executed, 
it  prompts  the  user  to  specify  the  input  geometry  file,  which  is  the  converted  facet  file  from 
atp.x.  If  no  input  file  is  available,  simple  geometries  such  as  quadrilaterals,  disks,  or  spheres 
can  be  generated  by  BUILDN5.X.  The  program  then  asks  for  a  set  of  simulation  parameters 
that  define  the  particular  calculations  to  be  performed.  At  this  point  the  excitation  voltages 
and  currents,  surface  impedances,  firequencies,  and  observation  points  are  entered.  The 
program  BUILDN5.X  conpiles  aU  of  this  information  and  constructs  an  ASCII  file  properly 
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formatted  for  use  in  PATCH.  This  file  must  be  renamed  inpatch  before  PATCH  is  executed. 
A  sample  inpatch  fiOie  for  the  mesh  in  Figure  2.8  is  listed  in  Appendix  A.  The  program 
BIJ1LDN5.X  is  also  capable  of  displaying  the  geometry  on  the  screen  or  writing  a  postscript 
file.  However,  geometry  display  requires  the  DISSPLA  graphics  software  package. 

When  PATCH  is  executed,  a  number  of  output  files  are  generated.  Outpatch  is  a  data 
summary  that  includes  a  listing  of  all  the  input  parameters,  geometry  (nodes  and  edges),  and 
the  computed  fields  and  currents.  Three  additional  files,  ampp.m,  ampt.m,  and  ang.m  are 
saved  with  .m  extensions  for  use  in  MATLAB.  These  files  contain  the  calculated  radiation 
pattern  data,  as  noted  in  Table  2.3.  MATLAB  is  the  primary  means  by  which  the  calculated 
pattern  data  is  displayed.  The  .m  files  are  simply  loaded  into  a  plotting  routine  to  view  the 
data  points. 

Other  PATCH  related  programs  include  bldmatx,  knit.x,  and  combine.in.  The 
program  bldmaLx  enables  a  display  of  the  mesh  model  in  MATLAB.  The  program  converts 
the  ASCII  facet  fiile  (from  atp.x  or  BUILDN5.X)  to  .m  files.  The  MATLAB  routine 
pltpatch.m  is  then  used  to  display  the  geometry.  Con[ibine.m  appends  one  geometry  file  to 
another  so  that  successive  nodes  and  edges  are  numbered  in  sequence. 

There  are  occasions  when  the  facet  files  generated  by  ACAD/AutoCAD  contain 
redundant  points  or  edges.  The  program  knitx  is  an  error  checking  routine  that  removes 
duplicate  edges  and  checks  for  other  file  errors.  This  program  provides  a  means  to  remove 
bad  data  and  generate  a  new  error-free  data  file. 
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D.  RADOME  PROGRAM 


The  FORTRAN  program  RADOME  is  used  to  compute  the  pattern  of  an  antenna 
radiating  in  the  presence  of  an  axially  symmetric  radome.  This  code  also  uses  the  method  of 
moments,  and  was  originally  developed  by  Francis  [Ref.  2]  and  later  improved  by  Klopp 
[Ref.  3].  See  also  [Ref.  7].  The  code  is  essentially  self  contained  and  does  not  require  the 
geometry  preprocessing  that  PATCH  does.  RADOME  is  restricted  to  axially  symmetric 
bodies  of  revolution  (BOR)  such  as  cones,  disks,  spheres,  hemisphares,  paraboloids,  cylinders, 
and  ogives.  These  basic  shapes  can  be  generated  internally  by  RADOME.  Other  generating 
ciuwes  can  be  specified  in  the  file  radinp,  which  is  simply  a  two-column  table  of  surface 
coordinates  p  and  z.  RADOME  produces  three  output  files  that  are  saved  with  a  .m  fiDie 
extension  for  plotting  in  MATLAB  (see  Table  2.3).  For  bookkeeping  purposes,  the 
underscore  character  in  the  filenames  is  occupied  by  a  letter  character  that  identifies  the  run. 
Just  as  in  the  case  of  PATCH,  each  observation  angle  in  ang_.m  has  a  corresponding  pattern 
value  in  eppol  .m  and  etpol_.m. 


Table  2.3  Output  files  for  PATCH  and  RADOME 


PATCH 

RADOME 

Description 

output  files 

output  files 

ang.m 

ang_.m 

Azimuth  observation  angles  (0  points) 

ampp.m 

eppol_.m 

E((.  (4)-polarized  E-field  component,  in  dB) 

ampt.m 

etpol_.m 

Eq  (0-polarized  E-field  component,  in  dB) 
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m.  COMPUTER  CODE  VALIDATION 


A.  APERTURE  MODEL 

In  order  to  validate  the  results  of  the  PATCH  and  RADOME  programs,  we  need  to 
accurately  model  the  antenna  used  in  the  anechoic  chamber  measurements.  Details  of  the 
chamber  facility  are  discussed  in  Chapter  IV.  For  all  of  the  experimental  data  that  was 
gathered,  a  Microline  56X1  Standard  Gain  Horn  was  used  as  a  receive  antenna.  This  horn 
has  an  e^qponential  taper  in  both  the  E-  and  H-planes  as  depicted  in  Figure  3.1.  Typically,  this 
type  of  horn  is  used  in  the  detection  and  measurement  of  microwave  power.  Additional 
specifications  for  the  horn  are  given  in  Table  3.1.  The  low  voltage  standing  wave  ratio 
(VSWR)  is  crucial  to  laboratory  measurements  since  it  prevents  undesirable  reflection  at  the 
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antenna  input  which  could  possibly  add  constructively  or  destractively  with  radome 
reflections,  thereby  producing  misleading  data. 


Table  3.1  Microline  56X1  Specifications 


Frequency 

Nominal 

Frequency  at 

Gain  Variation  over 

Accuracy 

Maximum  VSWR 

Range 

Gain 

Nominal  Gain 

Frequency  Band 

8.2 -12.4  GHz 

17  dB 

11.0  GHz 

+A  1.5  dB 

+/-0.3dB 

1.15 

1.  Patch  Code  Antenna  Model 

In  order  to  simulate  the  56X1  horn  in  PATCH,  a  model  of  the  antenna  aperture  was 
developed  based  on  slots  in  a  conducting  sheet.  The  excitation  amplitude  of  the  slots  is 
obtained  by  sampling  the  TEjo  mode  distribution  which  is  approximately  the  distribution  in 
the  horn  aperture.  A  rectangle  approximating  the  size  of  the  actual  horn  aperture  (5.4  x  7.3 
cm)  was  subdivided  into  a  12  x  13  grid  of  smaller  squares,  which  in  turn  were  divided  into 
two  faces  by  the  addition  of  a  diagonal.  At  10  GHz,  which  is  the  test  frequency  for  all 
measurements,  the  smaller  squares  are  0.2A  (0.6  cm)  in  dimension.  This  allows  the  individual 
element  to  be  electrically  small  while  keeping  the  grid  from  becoming  too  dense.  The 
aperture  grid  is  shown  in  Figure  3.2. 

The  next  step  is  to  determine  the  voltage  excitation  at  each  slot.  If  a  voltage  is  applied 
across  a  small  horizontal  slit  in  a  conductive  material,  the  resulting  far  field  pattern  will  have 
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Figure  3.2  Aperture  model  for  Microline  56X1  antenna 


a  vertical  E-plane  (y)  and  a  horizontal  H-plane  (x).  An  array  of  slits  with  varying  voltage 
amplitude  and  phase  excitations  can  be  used  to  control  the  array  far  field  pattern 
characteristics.  The  larger  the  array  becorres  in  the  horizontal  direction,  the  narrower  the  H- 
plane  beamwidth.  And  conversely,  the  larger  the  array  in  the  vertical  direction;  the  narrower 
the  E-plane  beamwidth.  For  the  Miaoline  antenna  model,  a  total  of  72  horizontal  edges  were 
chosen  as  slits.  These  are  denoted  by  small  dots  in  Figure  3.2.  The  extra  row  of  squares  at 
the  top  is  required  to  specify  the  polarization  of  the  slits  just  below.  The  size  of  the  aperture 
and  the  pattern  of  slits  were  adjusted  until  the  computed  pattern  agreed  with  the  measured 
56X1  pattern  in  gain  and  beamwidth. 

Because  of  TEjo  waveguide  excitation  imperfections  and  the  horn  flare,  the  aperture 
plane  of  the  56X1  does  not  have  uniform  illumination  with  constant  phase.  The  actual 
illumination  can  be  modeled  by  a  cosine  on  a  pedestal  function  as  shown  in  Figure  3.3.  The 
pedestal  value,  or  Ulumination  present  at  the  edge,  is  A,  and  the  width  of  the  aperture  is  a. 
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I(X) 


Figure  3.3  Cosine  on  a  pedestal 
Mathematically,  the  curve  is  expressed  by 


I{x)  =  A  +  B  cos  (3.1) 

with  A  +  jB  =  1.  The  illumination  tapering  in  the  56X1  is  a  result  of  the  horn’s  curved  shape. 
This  is  done  primarily  to  keep  the  input  VSWR  low  [Ref.  8]. 

An  aperture  phase  error  occurs  also  as  a  result  of  the  flare  geometry.  The  wave  fronts 
that  originate  at  the  waveguide  entrance  are  spherical  in  shape.  When  these  wave  fronts  cross 
the  flat  exit  plane  of  the  large  end  of  the  horn,  the  center  of  the  wave  is  ahead  of  the  outer 
edges.  Across  the  exit  plane,  the  phase  varies  as 


e 


(3.2) 


where  is  the  maximum  quadratic  phase  error  at  the  edge  of  the  aperture,  in  radians. 


24 


Together,  the  cosine  taper  and  phase  error  produce  an  illumination  in  the  horizontal  direction 
which  is  expressed  as 


I{x)  = 


A 


(3.3) 


Tapering  in  the  vertical  direction  is  expressed  by  substituting  y  for  x  and  b  for  a  in  Equations 
3. 1  through  3.3.  Since  the  taper  and  phase  error  in  the  vertical  direction  is  not  necessarily 
the  same  as  the  horizontal  direction,  the  constants  A  and  B  will  not  be  the  same,  and  ^  will 
be  replaced  by  i/fy . 

Applying  the  illumination  function  to  the  aperture  model.  Equation  3.3  is  sampled  at 
the  slit  locations  to  obtain  the  array  excitation.  Different  values  for  the  pedestal  and  phase 
error  were  investigated  until  there  was  satisfactory  agreement  with  the  measured  56X1 
antenna  pattern.  A  program  for  defining  a  two-dimensional  tapered  array  of  voltages  is  listed 
in  Appendix  B.  The  effect  of  increasing  the  taper  (B  >  A)  is  to  reduce  the  side  lobe  level  and 
increase  the  beamwidth.  The  effect  of  increasing  the  phase  error  is  to  fill  the 
pattern  nulls  between  the  main  beam  and  side  lobes.  Beamwidth  is  not  significantly  affected 
for  small  values  of  .  The  pattern  effects  are  illustrated  in  Figures  3.4  and  3.5. 

After  numerous  iterations,  the  best  agreement  with  measured  data  firom  the  56X1  horn 
occurred  with  E-plane  taper  values  of  A  =  0.6,  B  =  0.4,  and  =  35“  (denoted  by  0.6/0.4/35“). 
The  corresponding  best  H-plane  parameters  were  chosen  as  0.1/0.9/70“.  The  patterns  for 
both  planes  are  shown  in  Figures  3.6  and  3.7.  The  agreement  is  good  in  the  main  beam,  with 


25 


Effect  of  QuadraHe  Phase  Error  {E-Plane 


Figure  3.5  Effect  of  quadratic  phase  error 
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about  2  dB  difference  in  the  side  lobe  level  for  PATCH.  This  is  to  be  expected  due  to  the 
sanpling  nature  of  the  model  itself.  For  the  simulations  using  radomes,  the  aperture  model 
described  above  is  used  as  a  reference.  By  comparing  these  patterns  with  those  obtained  with 
the  radome  present,  the  effect  of  the  radome  can  be  deduced. 
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E-Plar«  Patt«rr  CempArlMn 


Figure  3.6  Comparison  of  main  beam  E-plane  patterns 


H’Pbjnd  Pattern  Comparison 


Figure  3.7  Comparison  of  main  beam  H-plane  patterns 
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2. 


Radome  Code  Antenna  Model 


The  code  RADOME  can  handle  circular  or  rectangular  antenna  apertures  of  arbitrary 
illumination.  The  antenna  pattern  is  computed  by  integrating  the  specified  distribution  at  a 
given  spherical  coordinate  location  (R,  0,  0).  No  approximations  are  made  in  the  calculation 
other  than  numerically  evaluating  the  integrals. 

When  using  a  rectangular  aperture,  the  RADOME  program  asks  for  the  antenna 
dimensions  in  wavelengths.  At  10  GHz,  the  56X1  horn  is  1.8  x  2.4  wavelengths.  In  order 
to  match  the  output  with  the  measured  data,  the  dimensions  and  distribution  were  adjusted. 
As  in  the  PATCH  model,  the  program  was  modified  to  accommodate  tapering  and  quadratic 
phase  error.  The  beamwidths  were  matched  in  RADOME  using  an  aperture  of  2.0  x  2.6 
wavelengths.  The  values  used  for  A,  B,  and  ^ ,  respectively,  are  1/0/50®  for  the  E-plane  and 
0.1/0.9/70®  for  the  H-plane.  The  resulting  patterns  are  shown  in  Figures  3.6  and  3.7.  There 
is  good  agreement  with  the  56X1  measured  data.  Furthermore,  for  the  purpose  of  radome 
evaluations,  there  is  satisfactory  agreement  between  the  patterns  produced  by  RADOME  and 
PATCH.  This  demonstrates  that  when  modeling  the  same  antenna  with  both  programs,  there 
is  some  correlation  in  the  required  tapers  and  phase  errors  and  similar  results  are  achieved. 
B.  SCATTERING  DISK  TEST  CASE 

The  data  in  Figures  3.6  and  3.7  show  that  both  the  RADOME  and  PATCH  antenna 
models  agree  with  measured  data  when  radiating  in  free  space.  When  a  radome  is  present, 
some  waves  scattered  by  the  radome  can  return  to  the  antenna  and  be  either  received  or 
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rescattered.  This  interaction  between  the  antenna  and  radome  can  affect  the  side  lobe  level 


of  the  transmitted  field,  especially  at  wide  angles.  The  RADOME  code  does  not  include  these 
multiple  reflections  between  the  antenna  and  radome,  but  PATCH  does.  The  strength  of  the 
antenna  reflection  in  PATCH  can  be  controlled  somewhat  by  varying  the  surface  impedance 
on  the  aperture  array  plate. 

A  good  test  case  for  evaluating  the  effects  of  multiple  reflections  is  to  place  a  perfectly 
conducting  (PEC)  disk  in  fi:ont  of  the  aperture  and  compare  the  results  for  both  codes.  The 
geometry  is  shown  in  Figure  3.8.  The  diameter  of  the  disk  is  2A,  and  it  is  positioned  5A  in 
front  of  the  aperture. 


Figure  3.8  PEC  disk  in  front  of  aperture 


The  resulting  patterns  are  shown  in  Figure  3.9.  The  patterns  from  the  main  beam  up 
to  the  fourth  side  lobe  agree  very  well  (0  <  80°).  The  PATCH  array  plate  impedance  is 

10,000  ohms,  which  makes  it  essentially  invisible  to  fields  scattered  from  the  disk.  Beyond 
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6=  90®,  PATCH  gives  a  number  of  additional  lobes  due  to  the  disk  scattered  field  adding  and 
canceling  with  the  array  field.  RADOME  has  no  antenna  field  in  the  rear  hemisphere,  and 
thus  the  pattern  in  this  region  is  determined  solely  by  the  disk  scattered  field.  Figure  3.9 
indicates  that  interactions  between  the  antenna  and  disk  are  not  significant  m  the  range 
0®  <  6  <  9Q  ,  which  is  expected  because  the  antenna  models  are  transparent  to  waves 
scattered  from  the  disk. 

In  order  to  test  the  multiple  reflection  effects,  PATCH  was  run  with  differentimpedance 
conditions  on  the  apertureplate.  Thesurfaceimpedanceoftheaperturewasvariedfrom  l,000to 
20,000 ohms,  with  and  without  the  diskpresent.  For  the  first  case  (aperture  alone),  the  results  are 
shown  in  Figure  3.10,  There  is  virtually  no  change  in  the  H-plane  pattern,  and  only  a  small  change 
in  the  E-plane  pattern  when  the  aperture  surface  impedance  is  varied. 

The  second  case  has  the  aperture  with  the  PEC  disk.  Again  the  surface  impedance  of  the 
aperture  plate  was  varied,  and  the  resulting  pattern  s  are  shown  in  Figure  3.11.  This  case  is  more 
interesting  in  that  the  surface  impedance  of  the  aperture  is  important  with  regard  to  multiple 
reflections  between  the  two  structures.  The  patterns  show  that  the  aperture  is  relatively 
transparent  for  impedances  of  1 ,000 to  20,000  ohms,  demonstrated  by  the  fact  that  these  curves 
are  grouped  together.  When  the  aperture  is  a  PEC,  we  see  the  pattern  diverges  from  the  others  as 
a  consequence  of  waves  being  reflected  back  and  forth  between  the  aperture  and  disk. 

In  summary,  the  PATCH  code  allows  a  wide  range  of  antenna  surface  impedances  which 
directly  affect  the  strength  of  the  interaction  between  the  antenna  and  scattering  body. 
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Figure  3.9  Patterns  with  aperture  and  disk 
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H- Plane  Pattern  (Aperture  Only) 


Figure  3.10  Effect  of  antenna  surface  impedance  on  the  aperture  pattern 
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Figure  3.11  Effect  of  antenna  surface  impedance  on  the  aperture  and  disk  pattern 
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IV.  DATA  ANALYSIS 


A.  ANECHOIC  CHAMBER  FACILITY 

The  Naval  Postgraduate  School’s  anechoic  chamber  was  used  for  all  measurements 
presented  in  this  thesis.  The  chamber  is  approximately  25  feet  long,  12  feet  wide,  and  10  feet 
high.  A  Hewlett-Packard  X890A  microwave  horn  serves  as  the  transmit  antenna  at  one  end 
of  the  chamber.  The  geometry  of  this  E-plane  sectoral  horn  is  shown  in  Figure  4.1.  A 
rotating  pedestal  is  located  approximately  19  feet  from  the  transmit  antenna.  The  pedestal 
holds  the  receive  antenna  (Microline  56X1)  and  test  radome,  and  is  capable  of  rotating 


Figure  4.1  HP  X890A  microwave  horn 

through  360“  in  azimuth.  Both  of  the  antennas  are  fed  by  coaxial  cable  via  WR-90  waveguide 
adaptors.  The  chamber  walls,  floor,  ceiling,  and  pedestal  are  covered  with  carbon-based  radar 
absorbent  material. 
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The  measurement  system  uses  a  miCTOwave  an5)lifier,  synthesized  sweeper,  frequency 
converter,  and  network  analyzer  to  generate,  detect,  and  process  microwave  signals.  The 
pedestal  servo  motor  and  network  analyzer  are  connected  to  a  personal  computer  that  serves 
as  the  main  controlleryinterface  for  the  chamber.  A  schematic  of  the  chamber  equipment  is 
shown  in  Figure  4.2. 

Measurements  were  taken  for  -90°  £6<  90° ,  with  ^=0°  pointing  directly  at  the 
transmit  antenna.  Data  was  not  recorded  for  angles  beyond  ±90° ,  since  the  receive 
antenna  faces  the  rear  of  the  chamber  in  this  sector.  In  practice,  the  radiation  pattern  will  be 
close  to  zero  because  of  the  missile  body.  Using  Labview  3.0  software  to  control  the 
instrumentation,  data  was  recorded  for  a  transmit  frequency  of  10  GHz  at  1°  increments.  The 
recorded  signal  is  an  average  of  32  samples  taken  at  each  azimuth  angle.  Sampling  reduces 
the  effect  of  noise  and  provides  more  repeatable  data. 

B.  PHASE  CENTER  FOR  MICROLINE  56X1  HORN 

In  the  far  field  of  a  microwave  horn  antenna,  it  is  desirable  to  assign  a  reference  point 
for  a  given  frequency  such  that  the  phase  of  the  radiated  wave  is  independent  of  6  and  (f>. 
This  point  is  called  the  phase  center  of  the  antenna.  When  viewed  from  a  far  field  observation 
point,  the  wave  fronts  radiated  by  the  antenna  are  spherical  waves  that  appear  to  originate 
from  the  location  of  the  phase  center  [Ref.  9]. 

A  graphical  method  of  obtaining  the  phase  centers  for  E-  and  H-plane  sectoral  horns 
is  given  in  [Ref.  9].  For  a  given  frequency  and  flare  angle,  the  phase  center  distance  is 
obtained  from  a  graph.  An  average  flare  angle  was  determined  for  the  E-  and 
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Figure  4.2  Anechoic  chamber  facility 
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H-planes  of  the  exponentially  tapered  56X1  horn.  Phase  centers  were  then  obtained  for  the 
horn  using  the  graph  in  [Ref.  9].  Using  these  estimates,  the  antenna  model  used  in  PATCH 
was  offset  from  the  physical  aperture  plane  (i.e.,  the  location  of  the  horn  aperture  when 
measurements  were  made)  by  the  phase  center  distance.  This  adjustment  is  shown  in 
Figure  4.3. 


The  E-  and  H-plane  phase  centers  are  an  approximation  since  the  56X1  horn  is  not 
a  true  sectoral  horn  in  either  plane.  To  obtain  the  best  phase  center  location,  the  PATCH 
antenna  model  was  stepped  through  different  locations  about  the  initial  estimate  until  the 
calculated  antenna  patterns  converged  to  the  measured  patterns.  This  location  was  then 
deemed  the  phase  center  for  a  given  plane.  The  final  values  for  the  56X1  antenna  are  0.0  cm 
for  the  E-plane  and  -3.0  cm  for  the  H-plane.  In  other  words,  the  E-plane  phase  center  is 
coincident  with  the  horn  aperture,  while  the  H-plane  phase  center  lies  3.0  cm  behind  it 
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C.  AGM-88  HARM  RADOME 

Pattern  measurements  were  made  on  an  AGM-88  High  Speed  Antiradiation  Missile 
(HARM)  radome  that  was  mounted  on  the  chamber  pedestal  in  front  of  the  Microline  56X1 
antenna.  This  set  up  is  shown  in  the  upper  part  of  Figure  4.4.  The  position  of  the  antenna’s 
aperture  plane  relative  to  the  radome  was  noted  so  that  an  adjustment  to  the  horn’s  phase 
center  could  be  made  in  the  computer  model. 

A  computer  model  for  the  HARM/56X1  antenna  configuration  mentioned  above  was 
generated,  and  is  shown  in  Figure  4.5.  The  antenna  model  was  located  at  the  E-  or  H-plane 
phase  center  for  each  respective  pattern  calculation. 


Figure  4.4  HARM  (upper)  and  AIM-9C 
(lower)  radomes  in  anechoic  chamber 
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Before  a  comparison  of  calculated  results  with  measured  data  can  be  made,  it  is 
important  to  examine  the  effects  of  both  radome  and  aperture  impedance  on  the  computer 
generated  results.  First,  the  effect  of  varying  aperture  impedance  was  examined.  This  is  the 
impedance  of  the  rectangular  aperture  plate  that  was  used  as  the  model  for  the  Microline 
56X1  antenna.  Figure  4.6  shows  the  E-  and  H-plane  patterns  for  aperture  impedances  of  100, 
500,  and  10,000  ohms.  The  E-plane  pattern  shows  the  most  change.  As  impedance  is 
inaeased,  the  side  lobe  peaks  tend  to  migrate  toward  the  center  ( ^  =  0®)  and  the  main  beam 
narrows  approximately  5°  at  -10  dB.  Side  lobe  levels  in  the  E-plane  are  decreased  with 
increasing  aperture  impedance.  The  H-plane  pattern  shows  relatively  little  change  with 
increasing  aperture  impedance.  The  main  beam  shows  virtually  no  change,  and  side  lobe 
levels  are  not  significantly  affected.  There  is  some  movement  toward  the  center  (as  in  the 


Figure  4.5  HARM  radome  model 
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E-plane)  in  the  second  and  third  side  lobes.  Thus,  the  value  of  aperture  impedance  determines 
how  transparent  the  antenna  is  to  waves  reflected  and  scattered  back  from  the  radome. 
Figure  4.6  shows  that  the  effect  of  aperture  irrpedance  on  the  HARM/56X1  model  is  minimal. 
There  are  no  gross  pattern  changes.  The  effects  are  mainly  in  the  low  side  lobe  levels  and  in 
azimuth. 

Secondly,  the  effect  on  radome  irrpedance  was  examined.  The  surface  impedance  for 
the  HARM  radome  is  278-731  ohms,  as  calculated  in  Chapter  in.  Pattern  calculations  were 
made  with  two  additional  values  of  radome  impedance,  208-/41  ohms  and  365-/54  ohms.  The 
effects  of  varying  this  impedance  are  shown  in  Figure  4.7.  In  both  E-  and  H-  planes  the  main 
beam  is  virtually  not  affected  and  there  is  a  slight  decrease  in  side  lobe  level  with  increasing 
radome  impedance.  This  indicates  that  the  more  transparent  the  radome  is  (i.e.,  the  greater 
the  impedance)  the  less  the  magnitude  of  reflections  between  the  antenna  and  radome. 
Overall  the  effect  of  radome  impedance  in  the  HARM/56X1  model  is  minimal,  with  only  low 
side  lobe  levels  being  affected. 

Now  that  the  effects  of  radome  and  aperture  impedance  have  been  examined,  a 
comparison  to  measured  data  can  be  made.  The  values  of  radome  and  aperture  impedance 
used  in  PATCH  were  278-/31  ohms  and  100  ohms,  respectively.  These  values  gave  the  best 
agreement  with  ensured  data.  The  patterns  are  shown  in  Figure  4.8.  The  E-plane  patterns 
show  the  same  general  agreement  in  side  lobe  structure  although  the  measured  side  lobes  are 
3-5  dB  higher.  The  main  beam  patterns  agree  well  down  to  -7  dB  where  the 
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Figure  4.6  Effect  of  varying  aperture  impedance 
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Figure  4.8  Comparison  of  measured  and  computed  patterns  for  HARM  radome 
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measured  beam  broadens  by  approximately  6®.  This  is  characteristic  of  a  large  quadratic 
phase  error  or  surface  wave  and  will  be  addressed  further  in  Chapter  V.  The  H-plane  pattern 
shows  good  agreement  in  the  side  lobe  structure,  and  better  agreement  than  the  E-plane 
pattern  in  the  main  beam  region. 

D,  AIM-9C  RADOME 

Pattern  measurements  were  also  made  on  an  AIM-9C  radome  mounted  on  the 
chamber  pedestal  in  front  of  the  Microline  56X1  antenna.  The  position  of  the  antenna’s 
aperture  plane  relative  to  the  radome  was  noted  so  that  an  adjustment  to  the  horn’s  phase 
center  could  be  made  in  the  computer  model. 

The  PATCH  computer  model  of  the  AIM-9C/56X1  antenna  configuration  is  shown 
in  Figure  4.9.  The  actual  radome  and  antenna  as  mounted  on  the  chamber  pedestal  is  shown 
in  the  lower  part  of  Figure  4.4.  The  antenna  model  was  located  at  the  E-  or  H-plane  phase 
center  for  each  respective  pattern  calculation. 

As  with  of  the  HARM  radome,  the  effects  of  aperture  and  radome  impedance  were 
examined  for  the  AIM-9C  radome  and  56X1  model.  Figure  4.10  shows  the  effect  of  varing 
aperture  irrpedance.  In  the  E-plane,  increasing  impedance  narrows  the  main  beam  and  alters 
the  side  lobe  pattern.  In  the  H-pIane  the  main  beam  is  largely  unchanged  while  the  side  lobe 
structure  is  markedly  changed.  This  was  a  key  factor  in  determining  the  best  choice  for 
aperture  impedance.  A  value  of  100  ohms  produced  the  best  agreement  with  measured  data. 

The  effect  of  radome  impedance  is  shown  in  Figure  4.11.  Both  E-  and  H-plane 
patterns  show  only  minimal  changes  as  the  impedance  is  varied  above  and  below  the  value  of 
234-yi66  ohms  that  was  determined  in  Chapter  III.  This  minimal  effect  is  most  likely  the 
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result  of  the  AIM-9C  radome  being  spherical  in  shape.  Most  of  the  reflected  waves  between 
the  aperture  and  radome  are  at  normal  incidence  to  the  radome. 

For  the  comparison  of  measured  and  calculated  patterns,  an  aperture  impedance  of 
100  ohms  was  used.  This  is  the  same  value  used  for  the  HARM  radome  case.  The  radome 
impedance  was  234-yi66  ohms.  Both  of  these  values  produced  results  that  matched  the 
characteristics  of  the  measured  patterns.  The  comparison  is  shown  in  Figure  4.12.  The  E- 
plane  patterns  have  the  same  basic  shape,  yet  the  measured  side  lobes  are  approximately  5  dB 
higher.  There  is  good  agreement  in  the  main  beams  of  the  H-plane  pattern,  but  some 
variation  between  the  side  lobe  patterns. 


Figure  4.9  AIM-9C  radome  model 
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Figure  4. 1 1  Effect  of  varying  radome  impedance 
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Figure  4.12  Comparison  of  measured  and  calculated  patterns  for  AIM-9C  radome 
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V.  CONCLUSION 


The  analysis  of  radomes  presented  in  this  thesis  uses  a  method  of  moments  (MM) 
solution  to  the  electric  field  integral  equation.  The  method  of  moments  provides  a  rigorous 
solution  that  is  well-suited  for  computer  coding.  The  MM  computer  codes  available  at  the 
Naval  Postgraduate  School  are  RADOME  and  PATCH.  The  primary  objective  of  this  study 
was  to  evaluate  PATCH,  since  it  provides  a  higher  degree  of  flexibility  and  accuracy  than 
RADOME.  The  collection  of  pre-  and  post-processing  programs  that  support  PATCH  is 
described  in  Chapter  H.  That  chapter  serves  as  a  condensed  “user’s  manual”  for  PATCH  that 
provides  sufficient  details  for  modeling  antenna  and  radome  problems. 

The  PATCH  code  antenna  model  was  validated  in  Chapter  HI  by  comparing  the 
computed  patterns  to  measured  data.  Once  validated,  PATCH  was  shown  to  be  more 
versatile  than  RADOME  in  a  number  of  areas.  First,  PATCH  can  model  any  arbitrary  radome 
shape,  while  RADOME  is  limited  to  bodies  of  revolution.  Second,  PATCH  takes  into 
account  all  multiple  reflections  between  the  radome  and  antenna.  The  calculation  of  this 
interaction  is  desirable  for  predicting  radome  performance  at  wide  angles  90°) . 

The  analysis  of  the  HARM  and  AIM-9C  radomes  demonstrate  the  versatility  of 
PATCH.  Models  were  developed  for  the  antenna  aperture  and  radomes,  and  a  range  of 
materials  (impedances)  were  examined.  Comparisons  of  the  calculated  patterns  with 
measured  ones  yielded  good  agreement.  The  pattern  differences  of  approximately  2-5  dB  are 
mainly  in  the  side  lobes.  The  results  of  the  analysis  with  PATCH  strongly  suggests  that 
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fundamentals  of  accurate  computer  modeling  exist.  However,  further  refinements  can  still 
be  made. 

There  are  several  possible  sources  of  error  in  both  the  computer  model  and  the 
anechoic  chamber  measurements.  Computer  model  errors  include:  (1)  approximate  surface 
impedance  of  the  radome,  (2)  inaccurate  modeling  of  the  radome  base  and  rim,  and  (3) 
inaccurate  antenna  modeling.  The  HARM  and  A1M-9C  radome  models  assume  an  average 
surface  impedance  over  the  entire  radome.  In  practice,  waves  from  the  antenna  will  have 
different  incidence  angles  along  the  surface  of  the  radome.  Since  surface  impedance  is  a 
function  of  incidence  angle,  the  impedance  should  vary  along  the  surface  as  well.  A  more 
accurate  model  of  the  radomes  would  include  this  impedance  variation  as  a  function  of 
position  on  the  radome.  It  is  expected  that  this  would  induce  focusing  of  the  main  beam  and 
alter  the  side  lobe  pattern. 

The  HARM  radome  has  a  hard  plastic  ring  which  was  not  included  in  the  antenna 
model  The  AIM-9C  radome  model  could  be  further  improved  by  including  the  metal  cylinder 
that  is  attached  to  the  dielectric  radome  shell.  This  cylinder,  or  ring,  is  bonded  to  the  base  of 
the  alumina  shell  and  serves  as  a  means  of  attachment  to  the  missile  body.  It  was  not 
removable,  and  thus  the  entire  assembly  was  moimted  in  the  chamber  (see  Figure  4.4). 
Although  the  cylinder  was  behind  the  antenna  aperture  and  out  of  the  field  of  view  out  to 
6=90° ,  back  scattering  from  this  piece  could  influence  pattern  measurements.  A  more 
accurate  computer  model  would  include  this  cylinder. 
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The  antenna  aperture  model  (Figure  3.2)  attempts  to  simulate  the  Microline  56X1 
antenna  as  closely  as  possible.  The  aperture  illumination  is  tapered,  with  quadratic  phase 
error  introduced,  and  the  aperture  plane  is  displaced  to  the  phase  center.  Despite  these 
adjustments,  some  minor  discrepancies  remain  when  comparing  the  antenna  pattern  radiated 
in  free  space  to  measured  data.  One  possible  source  of  the  differences  is  edge  diffraction. 

Insofar  as  measurements  are  concerned,  sources  of  error  include:  (1)  radome  mounts, 
and  (2)  imperfections  in  the  anechoic  chamber.  The  hardware  that  is  used  to  mount  the 
radomes  in  front  of  the  56X1  antenna  can  introduce  additional  scattering.  The  radomes  are 
supported  by  two  braces  made  of  phenolic  material.  Although  precautions  were  taken  to 
ensure  the  radomes  were  centered  in  front  of  the  antenna,  some  sagging  or  misalignment 
could  have  occurred,  which  results  in  pattern  asymmetries. 

Ideally,  an  anechoic  chamber  simulates  free  space  by  absorbing  all  waves  incident  on 
its  walls.  In  practice,  the  walls  will  scatter  and  reflect  some  of  the  incident  waves,  especially 
if  the  absorbent  material  has  been  damaged  or  is  deteriorated.  The  shape  of  the  chamber  also 
plays  as  a  role  in  measurement  error.  The  Naval  Postgraduate  School’s  chamber  has  as  a 
slanted  comer  wall  which  makes  it  asymmetrical  with  respect  to  the  chamber’s  centerline.  If 
radiation  were  to  be  scattered  from  this  wall,  pattern  asymmetries  would  arise.  The 
asymmetrical  nature  (with  respect  to  ^  =  0° )  of  the  measured  patterns  is  partially  due  to  the 
chamber’s  slanted  comer  wall. 

A  characteristic  of  composite  radomes  is  that  they  are  irdldly  anisotropic.  The 
dielectric  constant  is  polarization  dependent  because  of  the  sandwich  constmction.  For 
example,  the  HARM  radome,  which  is  a  fiberglass  composite  sandwich,  has  a  slightly 
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different  pattern  for  a  cut  measured  at  the  12  o’clock  position  than  it  would  at  the  3  o’clock 
position.  Although  the  HARM  is  as  a  symmetric  body  of  revolution,  very  small 
discontinuities  in  the  composite  material  can  lead  to  noticeable  pattern  variations  as  the 
radome  is  rotated  about  the  longitudinal  axis  [Ref.  10].  The  same  characteristics  have  been 
observed  on  the  AIM-9C  radome. 

Figures  5.1  and  5.2  show  the  effect  of  placing  a  radome  in  front  of  an  antenna 
aperture.  Figure  5.1  shows  the  effect  with  the  HARM  radome,  and  Figure  5.2  shows  the 
effect  with  the  AIM-9C  radome.  All  of  the  plots  are  calculated  patterns  from  PATCH.  In 
both  cases,  the  main  beam  is  largely  unchanged,  while  the  side  lobe  pattern  and  levels  are 
altered  significantly. 

This  thesis  has  established  that  PATCH  is  a  powerful  radome  analysis  tool.  PATCH 
is  highly  versatile  and  can  accommodate  any  radome  shape  or  configuration,  with  size  limited 
only  by  available  computer  memory.  Its  accuracy  is  mainly  limited  to  the  fidelity  of  the 
antenna  and  radome  models.  Future  work  should  concentrate  on  determining  the  level  of 
model  detail  that  is  required  for  a  specified  pattern  level  accuracy.  This  will  require  accurate 
surface  impedance  and  antenna  models.  In  the  case  of  antenna  scanning  the  surface 
impedance  will  change  because  the  angle  of  incidence  on  the  radome  changes.  Finally,  the 
gain  with  and  without  the  radome  should  be  computed.  In  general  this  requires  numerical 
integration  of  the  field  scattered  by  the  radome. 
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Figure  5.1  Effect  of  HARM  radome  on  antenna  pattern 
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Figure  5.2  Effect  of  AIM-9C  radome  on  antenna  pattern 
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APPENDIX  A.  PROGRAM  LISTING  FOR  IMPEDANCE  CALCULATION 


This  appendix  contains  a  listing  of  a  MATLAB  program  that  calculates  the 
transformed  free  space  impedance  through  a  three-layered  radome  wall  at  various  incidence 
angles.  The  real  and  imaginary  impedance  values  are  plotted  versus  incidence  angle. 

clear 

%  Programmers:  Prof.  D.C.  Jenn,  various. 

%  compute  equivalent  radome  surface  impedance 
%tincd=input {' enter  incidence  angle  (deg) ' ) ; 

%tinc=tincd*pi/180  ; 
f=10. ; 
start=0 . ; 
stop=:80 .  ; 
dth=l . ; 

N=round( ( stop- start ) /dth) +1 ; 

w=3e8/f/le9; 

for  m=l:N 

tincd=start+dth* (m-l) ; 
tinc=tincd*pi/180 . ; 

%  thickness  of  each  region 
dl= . 100* . 0254/cos (tine) ; 
d2= . 15* . 0254/cos ( tine) ; 
d3= . 05* . 02  54 /cos (tine) ; 

%  relative  dielectric  constants 
erl=3.5; 
er2=1.5; 
er3=:3.5; 
k=2*pi/w; 

%  assume  lossless  radome  materials  (alphal=alpha2=alpha3=0 ) .  Find  betas. 
al=:0.  ; 
a2=0 . ; 
a3=0. ; 

bl=k*sqrt (erl ) ; 
b2=k*sqrt (er2 ) ; 
b3=k*sqrt (er3 ) ; 
gl=al+j*bl; 
g2=a2+ j*b2 ; 
g3=a3+j*b3 ; 

%  characteristic  impedance  of  each  section 
z0=377 . /cos (tine) ; 

zl=z0*sqrt  (1 . /erl)  ;  %  assumes  theta__i  in  each  region 

z2=z0*sqrt (1 . /er2) ; 
z3  =z0*sqrt { 1 . /er3 ) ; 

%  impedances  looking  into  each  boundadry 
zc=z3* (z0+z3*tanh(g3*d3) ) / { z3+z0*tanh (g3*d3 ) ) ; 
zb=z2* (zc+z2*tanh(g2*d2) ) / ( z2+zc*tanh (g2*d2 ) ) ; 
za=zl* (zb+zl*tanh(gl*dl) ) / ( zl+zb*tanh (gl*dl ) ) ; 

%  reflection  coefficient  at  the  input 
%R(m)  =  ( (za-zO) / (za+zO)  )  ; 
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T  (m)  =tincd; 

Re  (m)  =abs  { za)  ; 

Iin(in)  =atan2  (iinag{za)  ,real  (za)  )  *180/pi; 
end 

subplot (211) ,plot (T, Re) , grid 
xlabel ( 'ANGLE  OF  INCIDENCE  (DEG)') 
ylabel ( 'MAGNITUDE  (OHMS) ' ) 
subplot (212) ,plot (T, Im) , grid 
xlabel ( 'ANGLE  OF  INCIDENCE  (DEG)') 
ylabel (' PHASE  (DEG)') 


APPENDIX  B.  SAMPLE  INPATCH  FILE  LISTING 


The  appendix  contains  a  listing  of  the  file  inpatch  for  the  quadrangle  shown  in 
Figure  2.8.  There  are  20  edges  and  43  nodes.  The  nodes  are  listed  first,  with  their  respective 
Cartesian  coordinates  (jc,  y,  z).  The  edge  segments  are  listed  next,  with  their  respective 
endpoint  nodes. 

The  last  block  of  data  defines  the  calculation  parameters.  For  this  case  there  are  two 
voltage  potentials  (between  nodes  4  and  5  and  nodes  7  and  8),  each  with  a  value  of  1.2  +  yO.9 
volts.  There  are  two  surface  impedance  faces  (1 1  and  12),  with  a  value  of  100  +  jlOO  ohms. 
Calculations  are  specified  for  Q°<d^9Q°  at  1°  increments.  The  observation  plane  is 
0  =  90°  and  the  frequency  10  Ghz. 


sample  inpatch  file,  comment  line 


20  43 


1 

0. OOOOOOE+00 

O.OOOOOOE+00 

O.OOOOOOE+00 

2 

0. OOOOOOE+00 

l.OOOOOOE+00 

O.OOOOOOE+00 

3 

O.OOOOOOE+00 

2. OOOOOOE+00 

O.OOOOOOE+00 

4 

0. OOOOOOE+00 

3. OOOOOOE+00 

0. OOOOOOE+00 

5 

1. OOOOOOE+00 

0. OOOOOOE+00 

0. OOOOOOE+00 

6 

l.OOOOOOE+00 

l.OOOOOOE+00 

O.OOOOOOE+00 

7 

1. OOOOOOE+00 

2. OOOOOOE+00 

0. OOOOOOE+00 

8 

1. OOOOOOE+00 

3. OOOOOOE+00 

O.OOOOOOE+00 

9 

2. OOOOOOE+00 

O.OOOOOOE+00 

0. OOOOOOE+00 

10 

2. OOOOOOE+00 

l.OOOOOOE+00 

O.OOOOOOE+00 

11 

2. OOOOOOE+00 

2. OOOOOOE+00 

0. OOOOOOE+00 

12 

2. OOOOOOE+00 

3. OOOOOOE+00 

0. OOOOOOE+00 

13 

3. OOOOOOE+00 

O.OOOOOOE+00 

O.OOOOOOE+00 

14 

3. OOOOOOE+00 

1. OOOOOOE+00 

O.OOOOOOE+00 

15 

3. OOOOOOE+00 

2.000000E+00 

O.OOOOOOE+00 

16 

3. OOOOOOE+00 

3. OOOOOOE+00 

O.OOOOOOE+00 

17 

4. OOOOOOE+00 

0. OOOOOOE+00 

0. OOOOOOE+00 

18 

4. OOOOOOE+00 

1. OOOOOOE+00 

0. OOOOOOE+00 

19 

4. OOOOOOE+00 

2.000000E+00 

0. OOOOOOE+00 

20 

4. OOOOOOE+00 

3.000000E+00 

O.OOOOOOE+00 
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K)  O  MOO 
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2  2  3 

3  3  4 

4  15 

5  16 

6  2  6 

7  2  7 

8  3  7 

9  3  8 

10  4  8 

11  5  6 

12  6  7 

13  7  8 

14  5  9 

15  5  10 

16  6  10 

17  6  11 

18  7  11 

19  7  12 

20  8  12 

21  9  10 

22  10  11 

23  11  12 

24  9  13 

25  9  14 

26  10  14 

27  10  15 

28  11  15 

29  11  16 

30  12  16 

31  13  14 

32  14  15 

33  15  16 

34  13  17 

35  13  18 

36  14  18 

37  14  19 

38  15  19 

39  15  20 

40  16  20 

41  17  18 

42  18  19 

43  19  20 

0  0 


4  5  1.20000E+00  0.90000E+00 

7  8  1.20000E+00  0.90000E+00 


11  l.OOOOOE+02 

12  l.OOOOOE+02 

false.  0  0 

1 

0.  90.0000  91 

0 


2.00000E+02 
2 . OOOOOE+02 


90.0000  90 


0000  1 
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false. 

0 

1. OOOOOE+IO 
-1.00000 
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APPENDIX  C.  PROGRAM  LISTING  FOR  VOLTAGE  TAPERING 


This  appendix  contains  a  listing  of  a  MATLAB  program  that  calculates  a  two- 
dimensional  array  of  voltage  excitations  across  a  rectangular  grid.  The  voltages  are  tapered, 
with  a  cosine  on  a  pedestal  distribution.  The  voltages  are  at  the  pedestal  value  at  the  edge 
of  the  grid,  and  highest  (1  volt)  at  the  center. 

The  variables  that  are  defined  at  the  beginning  of  the  program  specify  the  grid  size  and 
distance  between  individual  voltage  slits.  The  maximum  quadradic  phase  error  at  the  edge 
of  the  grid  is  specifed  by  px  and  py.  The  parameters  that  specify  the  cosine  on  a  pedestal 

taper  are  in  tapx  and  tapy. 

% 

%  Program  name:  TAPER. M 
%  Programmer:  LT  Scott  M.  Herzog,  USN 
%  January  1996 
% 

%  Program  for  indexing  segments,  nodes,  and  voltages 
%  for  use  in  INPATCH. 


clear 

a= 

0.078; 

%  width  of  waveguide  (meters) 

xl  = 

-0.030; 

%  x-coord.  of  center  of  first  segment 

dx= 

0.012; 

%  increment  between  segments 

xn= 

0.030; 

%  x-coord.  of  center  of  last  segment 

b= 

0.066; 

%  height  of  waveguide  (meters) 

yi= 

0.033; 

%  y-coord.  of  first  row 

dy= 

-0.006; 

%  increment  between  rows 

yn= 

-0.033; 

%  y-coord.  of  last  row 

segl  = 

42; 

Q. 

'h 

number  of  1st  segment  in  block 

dseg= 

2; 

% 

increment  between  segments 

segn= 

52; 

% 

number  of  last  segment  in  block 

jseg= 

30; 

% 

jump  in  segment  increment  (last  to  first) 

j  Iseg: 

=40; 

% 

jump  from  last  to  last  segment 

nodl  = 

2; 

% 

number  of  1st  node  in  block 

dnod= 

2; 

% 

increment  between  nodes 

nodn= 

12; 

'b 

number  of  last  node  in  block 

jnod= 

4; 

% 

jump  in  node  increment  (last  to  first) 

j  Inod: 

=14; 

% 

jump  from  last  to  last  node 

px= 

70; 

% 

phase  error  in  x-dir  (deg) ,  H-plane 

py= 

•n-n  . 

35; 

Q. 

'o 

phase  error  in  y-dir  (deg) ,  E-plane 

for  y=yl:dy:yn; 
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for  x=xl:dx:xn; 
n=n+l ; 

phx=exp{ j*px*pi/180* (2*x/a) ^2) ; 

phy=exp( j*py*pi/180* (2*y/b) ^2) ; 

tapx= ( 0 . 1+0 . 9*cos (pi*x/a) ) ;  %  H-plane  taper 

tapy=(0.6+0.4*cos(pi*y/b) ) 7  '  %  E-plane  taper 

v=tapx*tapy*phx*phy; 

volt (n) =v; 

end 

end 

n=0; 

for  y=yl:dy:yn; 

for  seg=:segl :  dseg :  segn; 
n=n+l; 

segment (n) =seg; 

end 

segl=segn+ j  seg; 
segn=segn+ jlseg; 

end 

n=:0; 

for  y=yl:dy:yn; 

for  nod=nodl :dnod:nodn; 
n=n+l; 

node (n) =nod; 

end 

nodi =nodn+j nod; 
nodn=nodn+ j Inod; 

end 

da  ta=:  [segment '  node'  real(volt)'  imag  ( volt )  '  ]  ; 
save  volt  n  data  -ascii 
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APPENDIX  D.  PROGRAM  LISTING  FOR  GEOMETRY  COMBINING 


This  appendix  contains  a  listing  of  a  MATLAB  program  to  combine  two  separate 
elements  of  geometry,  such  as  an  aperture  and  a  radome.  If  the  aperture  geometry  was 
created  separately  from  the  radome  geometry,  each  file’s  nodes  and  edges  are  numbered 
independently,  starting  with  the  number  one.  A  geometry  file  cannot  be  appended  to  another, 
until  the  second  file’s  nodes  and  edges  are  renumbered.  This  program  renumbers  the  nodes 
and  edges  of  the  second  geometry  file,  appends  it  to  the  first,  and  saves  the  combined  file 

under  the  name  duo. 

%  Program  name:  COMBINE. M 

%  Prograrnmer :  LT  Scott  M.  Herzog,  USN 

%  April  1996 
% 

%  Combines  two  independent  ASCII  data  files  from  the  output  of 
%  BUILDN5.X. 

% 

%  Re-indexes  the  node  and  edge  numbers  of  the  second  data  file 
%  and  appends  it  to  the  first. 

% 

%  Notes  on  usage: 

q, 

'o 

%  (1)  Delete  comment  line  in  both  input  files, 

%  (2)  Rename  first  input  file  FILE1.TXT,  typically  the  aperture. 

%  (3)  Rename  second  input  file  FILE2.TXT,  typically  the  radome. 

%  (4)  Output  file  saved  as  DUO,  without  a  comment  line, 

%  This  must  be  added! 

clear 

idl  =  f open ( ' f ilel . txt * , ' rt ‘ ) ; 
filel=  f scanf ( idl , ‘ %e ' ) ; 
id2  =  fopen ( ' f ile2 . txt ' , ' rt ' ) ; 
f ile2=  f scanf ( id2 , ’ %e ' ) ; 

nodesl  =  filel(l,l); 
edgesl  =  filel (2,1) ; 
nodes2  =  file2(l,l); 
edges2  =  file2{2,l); 

anl 
znl 
ael 
zel 
an2 


=3;  %  row  no.  of  first  element  in  node  block  in  filel. 

=  nodesl*4+2;  %  “  "  "  last 

=  znl+1;  %  row  no.  of  first  element  in  edge  block  in  filel. 

=  znl  +  edgesl'^3 ;  %  ”  "  ”  last  "  "  "  “  "  " 

=  3; 


%  read  in  first  file. 

%  read  in  second  file. 
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zn2  =  nodGs2*4+2; 

ae2  =  zn2+l; 

ze2  =  zn2+edges2*3 ; 

comb= [nodesl+nodes2  edgesl+edges2 ] ; 

Nl=reshape ( f ilel (anl : znl, 1) , 4, nodesl) ;  %  construct  node  block  matrix 

N1=N1 ■ ; 

El=rGshape { f ilel (ael : zel , 1 ) , 3 , edgesl ) ;  %  construct  edge  block  matrix 

E1=E1 ' ; 

N2=reshape ( f ile2 (an2:zn2,l) ,4,nodes2) ; 

N2=N2 • ; 


Nl(;,4)  =  Nl(:,4)  +  0.07699; 

N2(:,4)  =  N2(:,4)  -  0,1452; 
%N2(:,2:4)  =0.001*N2(:,2:4); 
%N2(:,2:4)  =  0 . 02 54*N2 { : , 2 : 4) ; 
%N2(:,2:4)  =  0 . 0205*N2 ( : , 2 : 4 ) ; 
N2(:,2:4)  =  0 . 0282369*N2 ( : , 2 : 4 ) ; 


%  translate  aperture  forward. 

%  subtract  offset  in  Z-axis. 

%  rescale  from  mm  to  m. 

%  rescale  from  in  to  m. 

%  rescale  from  'units'  to  m. 

%  rescale  to  m,  for  some  files 


E2=reshape ( f ile2 (ae2 ; ze2 ,1) , 3 , edges2) ; 

E2=E2 ■ ; 

N2  {  :  ,  1)  =:N2  (  :  ,  1> +nodesl;  %  changes  node  numbers  in  second  file. 

E2 ( : , 1) =E2 { : , 1 ) +edgesl ;  %  changes  edge  numbers  in  second  file. 

E2 ( : , 2 : 3 ) =E2 { : , 2 : 3 ) +nodesl ; 

save  duo  comb  N1  N2  El  E2  -ascii 
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